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Abstract: As a precursor to studying the bound states of multiple non-extremal black holes in 
A/" = 2, d = 4 supergravity, we investigate the stability of a near-extremal D0-D4 black hole in the 
probe limit, when the parameters of the black hole solution lie within a certain regime. We determine 
whether it is possible to form bound states of this "core" non-extremal black hole with BPS probe 
particles, and whether it is possible for the "core" black hole to decay by the emission of such BPS 
probes either to a local minimum of the probe potential, or spatial infinity. We first carry out a 
qualitative analysis of the probe potential to determine when quantum tunneling of probes from the 
^ ■ black hole is possible. We then find the wavefunction of the scattered probe by using the WKB 
d [ approximation to solve the Dirac equation in the black hole background, and use this solution to 
compute the tunneling amplitude. 
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1. Introduction 

A class of black holes that has not been extensively studied so far consists of the non-extremal 
black hole solutions to = 2, d = 4 supergravity[^, ^, |5), that are a generalization of the well- 
known extremal (both BPS and non-BPS) "attractor" solutions[^. These solutions consist of the 
background metric, together with a set of complex scalars and electromagnetic gauge fields, and 
can be completely characterized by the black hole's electric and magnetic charges, the parameter c 
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giving the deviation from extremality, and the value of the scalars at spatial infinity. In this work we 
study a particular class of these black hole solutions that have one electric charge Qq, one magnetic 
charge Pi, and one complex scalar field. Thus these solutions can be specified by four independent 
quantities: Qq, Pi, the extremality parameter c, and the value of the scalar field at spatial infinity. 
We package these quantities into two parameters, i and i-, that we use to classify the various 
regimes in which these solutions lie. 

The BPS attractor solutions are known to exist in multi-centered configurations that are sta- 
ble bound statesO, 13|. These multicentered configurations should remain valid, stable solutions 



even when deformed away from extremality[||, but fully backreacted, non-extremal, multicentered 
solutions have not yet been found. As the probe limit of BPS particles in a single-centered BPS 
background gives interesting insights into the fully backreacted supersymmetric multicentered solu- 
tions (such as the equilibrium distance between the centers) [0, 13, 12 1, it is reasonable to assume 
that studying the behavior of BPS probes in the background of a single-centered non-extremal black 
hole should provide clues to the existence and stability of multi-centered non-extremal black hole 
solutions. 

We thus consider a central "core" non-extremal black hole, that can be surrounded by BPS 
probes (in a multi-centered configuration with a large black hole at the center and the remaining 
black holes being small enough relative to the "core" black hole that they can be treated as BPS 
probes), and study the static potential of the probes in the black hole background. We assume that 
the probes are small enough that mutual interactions between them can be ignored. If the probe 
potential has a local minimum whose free energy is lower than that of the free energy at the black 
hole horizon, then the probe can form a stable bound state with the background black hole[^, |6| (and 
when fully back-reacted, this could give a bound state of two black holes.) A single black hole whose 
charge is equal to the sum of the "core" charge and the probe charge is unstable to the emission 
of such a probe to form this bound state, as the probe can tunnel through the potential barrier to 
the local minimum of the potential. If the probe has a lower free energy at spatial infinity than it 
does at the black hole horizon, then the "core" black hole is unstable, as it can emit probes that can 
tunnel through the potential barrier to escape to infinity. We cannot predict the endpoint of this 
evolution, as at some point the probe approximation (where the interaction between the probes is 
negligible) will become inapplicable: however, it is possible that the "core" black hole will continue 
to emit probes to infinity until we are left with a hot dilute gas. 

In this work we consider the emission of such charged probes from a certain class of non-extremal 
black holes in the near-extremal case, where the parameters describing the black hole solutions satisfy 
c ^ J- ~ We would like to know if these black holes are unstable to the emission of charged BPS 
probes, and if so, to determine whether these probes can form bound states with the background 
black hole, and to compute the tunneling amplitudes for this emission process. We first study the 
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qualitative features of the potential for a static charged probe in the black hole background, finding 
the maximum of the potential and the classical turning points. We also compute the rate for a 
charged particle to tunnel through the static potential between the classical turning points. This 
naive result for the tunneling rate through the potential barrier may also be applied to the case 
where the parameters of the black hole solution satisfy c <^ J- <gC J- . 

We then find the wavefunction for a charged probe particle in the black hole background, by 
solving the curved space Dirac equation for this background. We use the WKB approximation to 
solve for the radial part of the wavefunction. We then compute the amplitude for a charged probe 
particle to be emitted from the black hole via quantum tunneling through the potential barrier, by 
calculating the ratio of the conserved current density at spatial infinity and at the black hole horizon. 

This paper is organized as follows. In Section |2| we introduce the non-extremal black hole solu- 
tions that we will be studying, together with the necessary formalism for describing these solutions. 
In Section ^ we give the action for a charged probe particle in this black hole background, and carry 
out a qualitative analysis of the static probe potential. In Section ^ we set out the Dirac equation 
in the black hole background, and in Section |4.2| we carry out the full tunneling analysis by solving 
the Dirac equation. We conclude in Section |5|. 

2. Non-extremal black holes in A/^ = 2, d = 4 supergravity 

The action for the bosonic part of four-dimensional N = 2 supergravity coupled to massless vector 
multiplets takes the form: 



where the = 1, . . . , n) are the vector multiplet scalars, the (/ = 0, 1, . . . , n) are the vector 

field strengths, the Gi are the dual magnetic field strengths, and G^^^ = d^d^lC is derived from the 
Kahler potential 



where Qq is the holomophic 3-form on the Calabi-Yau manifold X. The normalized 3-form $7 = 



The lattice of electric and magnetic charges T is identified with {X, 12) , the lattice of integral 
harmonic 3- forms on X. In the standard symplectic basis, a charge V can be written as F = (P^, Q/), 
with magnetic charges and electric charges Qj. We can define a canonical, duality invariant, 
symplectic product (, ) on the space of charges, which is given by: 




(2.1) 




(2.2) 



(r,f) = p^Qj-Q,p^ 



(2.3) 
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in the standard symplectic basis. The moduh-dependent central charge Z{T, z) of F is given by: 

Z{T,z) = -e^l'^{T,Q) (2.4) 

All the coefficients of the Lagrangian can be derived from a single prepotential F{X), where the 
are projective coordinates such that = X^z^ and X^ is a gauge degree of freedom. In this paper 
we will consider prepotentials of the form: 

^W = ^^> (2-5) 

with D = 1, so that we have one scalar field, z, and two electromagnetic vector potentials, AP and 
(together with their duals, and -B^.) In this case the moduli-dependent central charge for a 
charge 7 = {P^,P^,Qi,Qo) is given explicitly by: 



2Y^D(Im2;)'^ \ o 2 J 

This restriction still allows us to consider a large class of black hole solutions, as it has been shown that 
the general case with an arbitrary number of n vector multiplets may be reduced to an effective theory 
with a single vector multiplet given by the prepotential (|2.5| ), by applying a suitable truncation]^. 

2.1 The Black Hole Solution 

Non-extremal black hole solutions to this theory were first found in |Q and Q, and further studied 
in and We consider the D0-D4 solutions of DO-charge Qo and D4-charge Pi, which may be 
described by a charge vector T = (P^,Qq). In analogy with the well-known extremal black hole 
solutions]^, a non-extremal solution can be given in terms of two functions Hq, Hi: 

Hq = sinh(cr -|- C2), Hi = ^- — — sinh(cr -|- C4) (2-7) 

c c 

where C2 and C4 are constants, c denotes the deviation from extremality, and r is an inverse radial 
coordinate such that r — t- 00 at the black hole horizon and r — )• at spatial infinity (in the extremal 
limit, the functions Hq,Hi are harmonic.) If Qq and Pi are both positive, then the solution is BPS 
in the extremal limit c — )• 0. If Qq and Pi have differing sign (where without loss of generality we 
can take Qq < 0), then the solution is non-BPS in the extremal limit. In this work we will restrict 
ourselves to the near-extremal regime where c ^ i- and c <^ i-. 
The black hole metric is given by[Q, ^]: 

ds^ = -e^Uir)^^2 ^ g-2C/(r) ( _^^dr'' + -^^dQl] (2.8) 

Vsinh CT sinh cr / 

where 

-2U 



e 



^HQHf. (2.9) 
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The scalar field z = iy is given by: 



y-r-§f, (2^10) 



and the electromagnetic vector potentials are given by: 



2 



5o = Qo(l-cose)#, Bi = -^Uc^ + ^-c\dt. (2.11) 



The integration constants C2 and C4 are given by: 

3/2 



where ?/o is the value of the scalar field at spatial infinity. Note that the solution is completely 

determined by the charges {P^, Qo), the extremality parameter c, and the value uq of the scalar field 

at spatial infinity. Thus four independent parameters are needed to specify the solution. In the rest 

of this paper we will refer to the parameters {P^,Qo), c, C2, and C4, but it should be kept in mind 

that one of these is redundant. 

We will find it convenient to define the parameters: 

sinh C2 , sinh C4 

C2 = , C4 = (2.13) 

c c 

The ADM mass M of the black hole can be read off from the metric: 



4' 



M = ^,/=^ + ^ + y=^ + % (2.14) 



where yo is the value of y at spatial infinity. The entropy of the black hole is: 

/ I 3/2 



and the Hawking temperature is: 



Th = ^ (2.16) 



In this work we will use the radial coordinate r defined by: 

r = (2.17) 
smh CT 

In these coordinates the metric becomes: 

ds^ = -e^Uir^f + e-2^('-) I , \. dr^ + r^d^l I (2.18) 
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And the functions determining the solutions are: 

1^1/ cosh Co L \ , ^ , / cosh ca L , ^ 

Ho = IQol i^-^ + 1 + ^ j , Hi = + 1 + ^ j (2.19) 

3. The Probe Action 

Before solving for the wavefunction of a BPS particle of charge 7 in the black hole background, it 
will be helpful to study the probe action for such a particle in this background, which is given byjl0[: 

S^ = -Jt,ds-^ y"(7,A^)dx^ (3.1) 

where is the moduli-dependent mass of the particle, given by: 

/x=|Z(7,z)| (3.2) 



The static probe action is: 

Sp = - J ^t-\J (7, At)dt (3.3) 



Vpdt 



where Vp is the static probe potential. This can be written in the form Vp = Vg + Vemi where 
Vg = e^\Z{-f, z)\ is the mass term, given byQ: 



and 14m is the electromagnetic coupling term, given by: 



We can read off the electromagnetic vector potential for this configuration from the full (i.e. non- 
static) probe action (|3.1|): 



4Qo VV ^0 J 4Pi Hf 
^, = i2£l(i _cos0) (3.6) 
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3.1 Qualitative Analysis of the Probe Potential 

In this section we carry out a quahtative analysis of the static probe potential in the near-extremal 
regime c <^ ^^^^ ^^"^ useful to define the following quantities: 



En 



1 go 3 

4Qo 4Pi 



3IQ0IIP1P 



go 

IQol 



+ 



1^(7,^)1; 



1 

16 



+ 



3pi 
3pi 



+ 



6|Qo 



Pi 



\Qo 



C2IQ0I C4IP1 



+ 



6|Qo|c2 



Po 

' \PI\ 



Po 



(3.7) 
(3.8) 
(3.9) 

(3.10) 



\Pl\C4: \C2\Q0\ C4IP1 

We will carry out a full scattering analysis only in the case where ~ Jj, which allows the calcula- 



C4- 

1 

<^2 



tions to be simplified. Note that in this case we have fiQ ~ /i^, and Eq 

The probe has zero potential energy at the horizon r = 0, while its potential energy at spatial 
infinity is 

1 go 3 pi 



Vr 



p\r=oo 



1^0 



AQ0C2 4PiC4 



+ cAi 



(3.11) 
(3.12) 



Tunneling of a probe particle to spatial infinity is only possible if this potential energy is non-positive, 

1. 1 

so that ficxi ^ Ai{Eq — c) (note that Eq ~ l/c2 >>> c, so that in order for tunneling to spatial 



infinity to be allowed, we must have Ai > 0.) We can see from (3.6) that the intersection product 
(7, r) between the probe and the black hole determines the magnetic part of the vector potential, 
not the electrostatic part: so the sign of the intersection product is independent of the sign of the 
static potential at spatial infinity. 

From (|]^-(|^ and ( p^ we see that if the black hole is BPS in the extremal limit, then 
> unless the extremality parameter c = and the intersection product (7, F) between 



V, 



p|r=oo 



the probe charge 7 and the black hole charge F is zero, in which case K 



p|r=oo 



0. In all other 



cases, the emission and absorption of probes is infinitely suppressed from black holes that are BPS 
in the extremal limit. Thus, apart from the special case of the emission of probe particles with 
(7, F) = from BPS black holes, we are only considering the scattering of probes from black holes 
that are non-BPS in the extremal limit. This result is very natural, as the fully backreacted, two- 
centered BPS black hole solution has an angular momentum proportional to (Fi,F2) where Fi and 
F2 are the charges of the two black holes, whereas a single-centered BPS black hole has zero angular 



momentum pi|, 13, |T^. Thus, if (7, F) 7^ 0, then emission of the probe 7 from the background black 
hole is forbidden by conservation of angular momentum in the extremal limit, when the background 
black hole is BPS. 
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We can study the shape of Vp in the two regions r <^ J- and r c. For r <^ J- we have: 



Vp\r<^J, = ^Mo - c^i ( ^ 1 + ^ - 1 1 (3.13) 



Solving ^\r=ro = gives: 



where 

(3.15) 

Since V"{ro) < 0, this is a maximum of the potential. And since ^oo < ^li^Q — c) when emission 
to infinity is allowed, and /x^ = /Iq, we have the following lower bound on (remembering that we 
are in the near-extremal regime CC2 1): 

«o > (3.16) 

~ CC2 (3.17) 
This also puts the following upper bound on tq: 

1 

ro<^<^^, (3.18) 
'-2 

so we see that the maximum of the potential is indeed in the region r <^ i in all the cases where 
emission to infinity is allowed. 
In the region r ;:g> c we have: 

Vp\r>^c= "^"l^^Z^'K cA, (3.19) 
1 + rc2 

so that 

which is < in all the cases where emission to spatial infinity is allowed. Thus we sec that with 
i ~ i-, the potential has one maximum at r = ro in the region r ~ c, then decreases continuously, 
tending towards the constant value Vp{oo) = ^^o — — c) at spatial infinity. The potential 

does not have a local minimum, and thus the "core" black hole cannot form a bound state with a 
BPS probe. A probe that is emitted will tunnel through the potential barrier and escape to spatial 
infinity. 
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The maximum value of the potential at rg is given by 



Emax = Vp{ro) = cAi{l - Oq) 



(3.21) 



Thus we can write the energy for a general scattered/emitted low energy particle as: 



e = cAi{l — (3ao) 



(3.22) 



where 1 < /? < The classical turning points for this particle are: 



c 



1 



(3.23) 



A low energy probe particle that has Poincare energy less than or equal to the maximum value fimax 
of the potential can be emitted from the horizon, tunnel through the barrier between the classical 
turning points r_ and r+, and escape to spatial infinity. Alternatively, a particle may come in 
from infinity and scatter off the potential, either being refiected back from the barrier, or tunneling 
through from the turning point r_(_ to r_ before falling into the horizon. If the background black hole 
is BPS, then the static probe potential Vp for a probe particle with (7,r) = is actually fiat for all 
r, so there is no barrier to emission or absorption. 

Note that Vp takes the same form as (|3.13D in the region r <^ J- when J- ^ J-. Thus, if ao < 1, 

* ' ■ L.2 '-'2 '^4 

then the probe potential in this case also has a maximum at ro with the same value emax- Although 
a complete scattering analysis for regions r ^ c is more complicated in this case, the tunneling 
rate between the turning points r± of this potential is the same as in the case J- ~ J- and may be 



calculated in the same manner, as we will see in Section 3.2 



3.2 The Tunneling Rate Through The Probe Potential Barrier 

Even without completing the full scattering calculation, we can obtain an estimate of the tunneling 
rate Ft through the potential barrier, which is defined as: 



-/;+lp;idr' 



(3.24) 



where r± are the classical turning points for the probe in the static potential, and pr is the radial 
canonical momentum, given by: 



or 

The Lagrangian density £. can be found from the probe action (|3.l| ): 



(3.25) 



\ 



o2U 



-2U 



1 + 



At, 



(3.26) 
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where At is given by (|3.6| ). This gives 



fie ^^r 



+ -e-2^fl + 



(3.27) 



The conserved canonical energy e satisfies: 



e = - £ (3.28) 



C2 



+ W ( 1 + ^ ) e4 V + e^^fi^ (3.29) 



Thus we find that 



2 ^ 



1 + ^ 



((6 - Atf - e'^fi') . (3.30) 



For the range of parameters that we are considering for the background black hole, the classical 
turning points r± for the scattered probe particle always lie in the region r <^ c. Thus, defining the 



coordinate z = r /c , the classical turning points lie in the region z ^ 1 and are given by ( 3.23 ): 



z± = g(±l + ^(l-a2)(l-i^)) -1 (3.31) 

4/3 



We find that 



and so 



z+-z_ = ^V(l-a^)(/32-l) (3.32) 



lim |Re {i [pr{z+ - () - Pr{z- - C)]) I = -aiQovr(/3 - 1), (3.33) 



Tt = e-'^i°o^(/^-i) (3.34) 

As mentioned at the end of Section |3.1| , Fj also gives the tunneling rate through the potential barrier 
in the regime J- <^ as the probe potential V„ has the same form in the region r <^ J- as when 



4. The Dirac Equation 

In this section we carry out a full analysis of the scattering and emission of a charged probe of charge 
7 from the black hole. We would like to solve for the wavefunction of a probe of Poincare energy 
e scattering off the background black hole. In Q, it was found that the ground state wavefunction 
describing a light BPS probe particle of charge 7 in the background of another, heavy BPS particle 
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of charge F, is given by a monopole spherical harmonic|14| corresponding to a configuration with 
total angular momentum ((7,r) — l)/2. This can be thought of as the light BPS probe going into a 
spin-1/2 state aligned with the radial magnetic field of the background BPS particle, thus minimizing 
the energy of the configuration, and contributing one spin quantum opposite to the intrinsic field 
angular momentum, which is (7, T) /2. This is analogous to the problem of scattering a Dirac particle 
of charge Ze in the background of a magnetic monopole of strength g^^, where the intersection 
product (7, r) corresponds to the quantity Zeg. Since we are trying to solve the curved space 
version of this problem (though for higher energy states as well as the ground state), it is reasonable 
to assume that the wavefunction describing a probe particle emitted from the background black hole 
will obey the curved-space Dirac equation in the black hole background, and correspond to the probe 
going into a spin-1/2 state aligned with the radial magnetic field of the background black hole. 
Using the vierbein formalism, the Dirac equation in curved space is: 



(4.1) 



where /i is the mass of the particle given by ;U = \Z('y, z)\, is the electromagnetic gauge potential, 
and ^' is a 4-component spinor. In our case is given by ( |3.6D . Note that any choice of the vector 
potential A^p around a magnetic monopole must have singularities. Thus, the specific form of A(f, 



given in (3.6) represents a gauge choice that is non-singular in some region Ra (in this case, the 
region r > 0, < < 2-it, and < 9 < vr.) In order to cover the entire space outside the magnetic 
monopole, we must divide the space into two regions and use a different gauge in each region. We 
can define the second region as r > 0, < < 27r, and < 6* < vr, and use the gauge choice 

(7,r) 



A, 



-(1 + COS I 



(4.2) 



in this region. 

Because of this fact, the components of the Dirac spinor 



(4.3) 



in a magnetic monopole background are given by sections on a line bundle, not a function|14|. The 
angular part of the section can be expanded in monopole spherical harmonics Yj i „^{9, (j)) character- 
ized by the quantum numbers {j,l,m), where 



. = o,ii,... 

^ = |j|,|j|+l,|j|+2, 
m = —I, + 



(4.4) 

(4.5) 
(4.6) 
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and 



Yj,i,rn{e^ 0) = e^«'"+^)'^e,, i,„,(0) in region Ra (4.7) 
y,- 0) = e'"'"-^'^'^e,- in region R,, (4.8) 

for the same function Qj^i^mid)- 

For a suitable choice of vierbeins (for details, see Appendix we substitute the ansatz 



V'l = Ri{r)Yg-i/2,iA(^, 4>)e"^^e-'^' (4.10) 
^^2 = R2{r)Y^+i/2,iA(^, 0)e^^*e--* (4.11) 
V's = i?3(r) Vl/2,^,„^(e, 0)e^^*e--* (4.12) 



V'o = iio(r) Vl/2,^,„^(e, 0)e^^*e--* (4.9) 

^q+l/2,l,m{ 
^q-l/2,l,m{ 

into ( [4.ip , where we have defined 

(4.13) 

In the case when g > 0, there is a possible solution with V'o = "02 = 0, and / = q — \- When 
g < 0, there is a possible solution with ■0i = ^3 = 0, and / = — (7 — ^. As expected, these solutions 
correspond to the probe particle being aligned with the radial magnetic field of the background black 
hole. The two cases are exactly analogous so from now on we will assume that (7 < 0. We obtain 
coupled radial equations for -Ro('') find R2{r) of the form: 

3U 

e'^idt + iAt + ifi)Ro + —drire^'^ R-?) = (4.14) 

r 

e-^{dt + iAt - in)R2 + —dr{re'^Ro) = 0. (4.15) 

r 

4.1 The WKB Approximation 

We cannot solve the radial part of the Dirac equation exactly in this background. So in order to 
compute the tunneling amplitude, we will solve for the radial components of the wavefunction using 
the WKB approximation. 

In order to apply the WKB approximation, we first substitute q ^ q/h and e — )• e/h. Using 
the ansatz (|4.9D and defining Tq^2{i^) = Te~^/^i?o,2(?') gives the following equation for To(r) (with an 
analogous equation for T2(r): 



° (ee-u - Ate-U + fi) dr 



Vl + c2/r2 



(4.17) 
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where a prime denotes differentiation with respect to r. Substituting the ansatz Tq^2 = -60,26*'^°'^''^ 
and taking terms to leading order in 1/h, we find: 

= ((« - - ■ (1-18) 

Note that this is the same expression that gives the canonical momentum in ( |3.30| ) . We need to go 
to the next order in 1/h to find the equation for BQ{r) (an exactly analogous equation gives B2{r)): 



BqSq + 2BqSq 



+ e^Vl + c2/r2- 



{ee-u - Ate-u + /i) ^ ' dr \ + 



BoS'o (4.19) 



This next order in the WKB approximation is necessary in order to derive the connection formulae 
needed to extend the solution past the classical turning points, where the WKB approximation 
becomes invalid. The connection formulae are derived in Appendix |y. We can solve equations 
( [4.18D -( p:.19| ) in the region r <^ i- (recall that both turning points lie in this region for the class of 
background black holes we are considering), and then patch this solution to the WKB solution in 
the region c <^ r. The details of the calculation are given in Appendix |^. 

4.2 The Tunneling Amplitude 

In the regime c <^ J- ~ J- we find that 

cosh C2 , cosh C4 1 (4.20) 
sinh C2 f« CC2 1 (4-21) 

In order to calculate the tunneling amplitude for a probe particle, we want to find the equivalent of 
the probability density of the particle wavefunction at different values of r. In the case of a Dirac 
spinor, this is given by the time component of the conserved current density J^, which is given by: 

= ^Y^^,^ (4.22) 
V^J'^ = 0, (4.23) 



where 



r := KTt", (4.24) 
^ := ^'"^7°, (4.25) 



and 7° indicates the flat-space gamma matrix. 

For the background metric ( 2.1^ ), is given by: 



j" = e-^(|i?oP + |i?2|'). (4.26) 
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We now want to compute the ratio of charge densities at spatial infinity r — )• c«, and at the black 
hole horizon r — )• 0. Given a conserved current J'^ satisfying V^J^ = 0, we can define the one- form 
= g^uJ'^ and write the conservation condition as: 

d(*J) = (4.27) 

We then define the charge passing through a hypersurface T-L via: 

Qn = - [ (4.28) 
Jn 

We take 7^ to be a hypersurface of constant time, t. We can then write: 



Qn = - [ d^xv^n^J/^ (4.29) 
Jn 

where hij is the spatial metric and = gfj,un'^ where n'^ is the normal vector to the hypersurface. 
We therefore have 

n*" = (1,0,0,0) (4.30) 
n;,= (-e^^, 0,0,0) (4.31) 



p-GU 

\h\ = J . , , . .. r^sin^g. (4.32) 



(l + c2/r2) 



We can then write: 



Q-H= [ d^xy^\e^^J^ (4.33) 
Jh 

r I P-6U 

= / drddd^e^i ^-^r^sin^e (|i?op + |i?2p) (4.34) 

Jn V (1 + c7r2) 



2U 2 

drded<P^^Z^ {\Rof + \R2f). (4.35) 
H y'l + c^/r^ 



It follows that to get the charge density q{r) at r, we should integrate just over 9, 



q{r)= [ dg^ ^'^^'fr! i\Ro\' + \R2\') (4.36) 
JH Y^l + c2/r2 

~-^==(|i?oP + |i?2p) (4.37) 

We can now calculate the charge densities at the black hole horizon, and at spatial infinity, using 
the form of the Dirac spinor found in Appendix ^. We impose boundary conditions so that the 
wavefunction is entirely outgoing at r — )• oo. In order to obtain the emission amplitude, we should 
compute the charge density at r — t- using only the outgoing component at the black hole horizon. 
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We find that: 



(?(r)™ = r'{\Ro\' + \R2n = \Bo\' (4.38) 

Eor,,-r, ,2 , ID |2^ _ "0 |6 |2 /^T. 1 



for some constant Bq, where the tunnehng rate Tt is given by (PI). We can obtain the tunnehng 
amphtude by taking the ratio of the charge densities: 

= ^^(l-/^«o) I (4.41) 

ao (g-aiao7i-(/3-l) _ 1 gaiao7i-(/3-l)) ^ 

We can see that the tunnehng amphtude depends on the parameter ao, which is related to the 
difference between the gravitational and electrostatic parts of the static probe potential, as well as 
the energy of the probe (which is related to (1 — /3ao)). It also appears that the tunneling amplitude 
decreases to zero as the background black hole becomes extremal. 



5. Conclusion 

In this paper we have initiated a study of the stability of a class of non-extremal black holes in 
A/" = 2, d = 4 supergravity. Our results provide some clues as to how such a black hole perturbed 
away from extremality might evolve over time. If the black hole is BPS in the extremal limit, then 
for the range of parameters we have considered, then it can only emit probe particles 7 such that the 
intersection product (7, T) between the probe charge and the black hole charge is zero. The static 
potential felt by such a probe is flat, so that there is no barrier to emission or absorption. Thus, we 
would expect that such a black hole would decay by the emission of these probes to spatial infinity 
until it became extremal. 

If the black hole is non-BPS in the extremal limit, then it can emit probes of different charges. 
We found that the decay rate increases as the parameter Uq decreases, where is related to the 
difference between the gravitational and electrostatic parts of the static probe potential. Thus, we 
expect probes with small to be emitted preferentially by the black hole. For the class of black 
holes we studied, we found that the probe potential has a single maximum at r = ro, and decreases 
steadily as r — t- 00. Thus, as there is no local minimum in the probe potential, we expect that any 
particles emitted from the black hole will simply be ejected out to spatial infinity. 

There are many possible routes for further investigation: it would be useful to extend these 
results to arbitrary values of the parameters c, i, and i-, so that we are no longer restricted to 
the near-extremal limit. Another natural extension of our results would be to study the scattering 
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of fermionic particles in the black hole background, and to compare the relative emission rates of 
different types of particles from the black hole. 
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A. The Dirac Equation in the Vierbein Formahsm 



Here we give the explicit form of the Dirac equation for a particular choice of vierbein. We use the 
following representation for the 7-matrices in the Dirac equation (4.1): 



7 



1 
-1 



7 



-a' 
a' 



(A.l) 



and the following representation for S"'': 



(A.2) 



We use the following choice of vierbein: 



V3 



e ^rsinl 



1 + ^dr 



(A.3) 
(A.4) 
(A.5) 

(A.6) 



The Dirac equation for a spinor 



( V^o^ 



(A.7) 
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then becomes: 



3t/ , „ 



3U 



e"^ {dt + iAt + - ^5^(re-f V3)\/ 1 + ^ 



3t/ 



e-^(5t + iAt - i^)'^2 + —dr{re-^ipo)\ 1 + 



c2 



3U 



e~^idt + iAt - iiJ,)ip3 —dr{re~2'ilj^)Jl + — 



where we have defined 
We now use the fact that 

V = -de-^d^+jcot9 

smt) 

_ 2 

V = -de + ^—^dA> - j cot 6 

smU 

are raising and lowering operators for the monopole spherical harmonics that satisfy: 



'DYj,l,m = -^/il+m-j + l)YJ_,,l,r, 



We substitute the ansatz: 



V'i=i?i(r)y,_i/2,i,^(0,0)e^«-^e--* 
^2 = R2{r)Yg+,/2,i,m{e,<f>)e'"^e-''' 
V'a = i?3(r)y,-i/2,/,™(^,0)e'^'^e--*, 

and see that ii q > 0, then there is a possible set of solutions to the Dirac equation with ■i/'o = V''2 = 0, 
and I = q — ^. Similarly, if g < 0, then there is a possible set of solutions with V'l = V's = 0) 
l = -q-i 
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B. Solving the Dirac Equation Using the WKB Approximation 



Here we give the details of the calculations for solving the Dirac Equation using the WKB approx- 
imation. As shown in Section 4.1, this comes down to solving equations ( 4.18| ) and ( 4.1S| ) for the 
radial parts of the Dirac spinor components. 
In the region r ^ c, Eq. ( [4.18 ) becomes: 

2 



S'i = EoAjal 1 + ^ 



(B.l) 



and Eq. ( [4.19| ) becomes: 



r 1 + 



E, 



1/2 



These equations can be solved exactly to give: 



(B.2) 



As r — )• oo, the solutions become: 



Bo{r) = Bo 

5o(r) = ±^(r + 4logr) 



——Bo{r) = — 



aiao 
S±{r) = ± — —r 



(B.3) 
(B.4) 



(B.5) 
(B.6) 



Imposing outgoing boundary conditions at infinity (so we only keep the far region solution is: 

r 

In order to patch the solutions in the regions r c and r <^ we evaluate Bo{r) and S-^-{r) in 
the limit r <^ i-, which gives: 



—-Bo{r) = 

r2Ei 
5+(r) = aiaologr 



(B.8) 
(B.9) 



Thus the solution in the region r c is: 



Ae^ 



(B.IO) 
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Defining z = in the region r ^ Eq. ( [4.18| ) becomes: 



V2 



4z2(l + 



((/3qo - ^YT~zf + (ag - l)z) , 



(B.ll) 



where a prime denotes differentiation with respect to z. And Eq. ( 4.19| ) becomes: 



2(l + z)5o(5o + 2z5o) + 2z l-2(l + z) 

+ 8z(l + z)So'5o = 
Equation ( |B.11| ) can be solved exactly to give: 
5o(z) = 5±(z) 



-Bo 5*0 



(B.12) 
(B.13) 

(B.14) 



ai 



-2/3aoArctanh(Vl + ^) - log(-2;) + 2ao log [-/3 + aoVl + z + /(z)] 
+ (1 - /3ao) log [1 - /(z) + ao (ao(-l + VTT^) + - + - /(z)))] 

+ (1 + /3ao) log [1 + /(z) + ao (-ao(l + ^^TT^) + /3'ao + /?(! - ^f^z + /(^)))] 

(B.15) 



where 



/(z) = V 1 - 2/3ao\/rT^ + (z + /32)a2^ 



(B.16) 



The form of Bq(z) may also be determined exactly, but we will not give it here as it is extremely 
complicated and not particularly illuminating. 

In order to patch the solution to the region r ^ c, we solve Eq. (p^-(Pl) for z 1: 



2^,2 



1"0 



Q/2 

^0 "1? 



4zSo('5o + zS'f^) + Sz'B'qS'q = 



to give: 



-Bn(z) 



S±{z) 



1 

Z4 



aioo , 
± log z 



Matching to ( [B.10|) , we find that the solution in the region z ^ 1 is: 



.laiao 



^4 1 1 

E^c^z^ 



z*^ 



(B.17) 
(B.18) 

(B.19) 
(B.20) 

(B.21) 
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We can extend this solution up to the outer turning point at z = z+, hy taking S+{z) to select the 
wave that gives this form of -0 for z 3§> 1. At the turning point, the WKB approximation becomes 
invalid, as Bq(z) diverges. Thus, in order to continue the WKB solution past 2:+, we must derive 
connection formulas that allow us to patch together the solutions in the regions z < z_ and z > Z-. 
We then repeat the procedure at the other turning point This can be done using Airy functions: 
detailed calculations are given in Appendix ^ 

Using the connection formulas, we find that the solution in the region z < z_ is: 

where the tunneling amplitude Tt is given by ( |3.34| ). The functions Bo{z) and S±{z) in the region 
z <^ 1 are determined by the equations: 

S'o' = ^{1 - Paof (B.23) 

450(5*0 + zS;;) + SzB'^S'o = 0, (B.24) 

which give: 

1 

2 Jaiao 

Bo{z) = ^^Bo (B.25) 

(1 -/3ao)2 E^C2 

So{z)±=±^{l-(3ao)logz (B.26) 

where Bq is a constant. The wavefunction in the region z <^ 1 is thus: 

pU/2 

^ ^ 1 Boe'^o/^ (B.27) 



r 



a. 







,^5o.-V4 ((r, - ^) .+^(1-^-) - . (r, + ^) .-^a-/^-)) (b.28) 



(l-/3ao)2 E^c 

It may appear that the solution blows up and becomes unphysical as z — )• 0, but as we will see in 
Section ( [4. 21) , this is not a problem as the charge density always remains finite. This solution is the 
limit of the general solution for z <^ 1. 

C. Deriving the connection formulas for the WKB approximation using Airy func- 
tions 

Here we outline the calculation for deriving the connection formulas that allow the WKB solutions 
to be extended across the turning points at 2:+ and z_. The general procedure is to solve radial 
equation using a linear approximation to the potential at the turning point, and then match this 
solution to the WKB approximation to the left and right of the turning point. 
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If we define the potential 



(/3ao - VTTzf + {al - 1)^ 



4z2(l + z) L 

then the radial equations for the WKB ansatz ip = Br{r)eh^^^'^^^''^^Yg^i^m{d , (I)) 
S'^ = V{z) 



2{l + z)Br{S' + 2zS") + 2z l-2(l + z) 
+ 8z{l + z)B'^S' = 
We can rewrite the second equation as: 



an 



BrS' 



2(1 + z)Br{S' + 2zS") + 8z(l + z)B'^S' + 2z[l- 2(1 + z)Miiz, ao,^)] BrS' = 
for some function Mi{z,ao, /3). 



(C.l) 

(C.2) 
(C.3) 
(C.4) 

(C.5) 



The exact radial equation, given by ( 4.16 ), is: 

T" + M{z, ao,P)T' + V{z)T = 



(C.6) 



for some function M{z,ao, (3), where T{r) = re ^^'^R{r) and R{r) is the radial part of the Dirac 
spinor components. 

C.l The Linear Approximation 

Consider a turning point at z = zq. Defining the coordinate w = z — zq, we can take the linear 
approximation to the potential around zq: 



Viz)^Vizo) + V'izo)w 
= pw 



(C.7) 
(C.8) 



where p = V'{zo). Note that we have /) < for = Z- and p > for = The radial equation 
can then be approximated as: 



T"{w) + M{zQ, ao,P)T'{w) + pwT{w) = 
This equation has the general solution: 

— wM(za .aci ,3) ^ 

T{w) = e [CiA\{w') + C2Bi(u;')] 

where K\[w') and Yii{w') are Airy functions, and we have defined 

M(zo,ao,/3)2 



w = {—p) + 



A(-p)i 



(C.9) 



(C.IO) 



(C.ll) 
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At each turning point we assume that there is a region in which the hnear approximation is vahd, 
and that the second term in the expression for w' is neghgible compared to the first, so that we can 
take: 



w' ^ {—p)3w. 



(C.12) 



C.2 The turning point at z+ 

First consider the turning point at z = z^. The general solution in the region around z^ where the 
linear approximation is valid, is: 

— wM (z^ ,aQ,,3) 



CiAi{w') + C2Bi{w'] 
where w' is defined as in ( p.l2| ). Since p is positive, we find that 

/ / N i i 

W ~ {—l)3p3W, 



(C.13) 



(C.14) 



and thus in order to match the above solutions to the WKB approximation to the left and right of 
2+, we need the asymptotic behavior of the Airy functions for w' ^ {—1)3 x — oo: 

1 



Ai(u'') 
Bi(u;') 



1 



(C.15) 
(C.16) 



and w' — )■ (—1)3 X oo: 



Ai(u;') 



1 2^'3/2 

-re3^ 



1 

r- /i 



^ 2 ,3/2 -2^/3/2 

-63 +e 3"^ 
2 



(C.17) 
(C.18) 



We now need to solve the WKB equations to the left and right of the turning point. The radial 
equations for the WKB ansatz in the region where the linear approximation is valid are: 

S'(wf = pw (C.19) 
2(1 + z+)Br{S' + 2z+S") + 8z+{l + z+)B'^S' + 2z+ [1 - 2(1 + z+)Mi(z+, ao, /?)] BrS' = (C.20) 



These equations have the solutions: 



and 



2i 



S{w) = ± — y^pw^ 
o 



Br 



BqC 2 



(C.21) 



(C.22) 
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where Bq is a constant and M(z-^-,a, (3) is the same function as in (103^) . Note that pw^ is negative 
on the left of the turning point, and positive on the right. 

Since we are imposing the boundary condition that only the outgoing wave is present at spatial 
infinity, the WKB solution to the right of the turning point is: 



—e 3 



1 



' fyW 



In order to match the linear solution in ( C.13| ) to the WKB solution, note that: 



2. ,n3 



2i , o 1 



(C.23) 



(C.24) 



and use the asymptotic form of the Airy functions as w' — )• (— 1)^ x oo. We assume that there is a 
region in which the linear approximation is valid, and \w'\ is large enough for the linear solution to 
be well approximated by this asymptotic form. We find that: 



Ci = 2/^(-l)T2pT2 5o 
C2 = 0. 



(C.25) 
(C.26) 



We now want to match the linear solution to the WKB solution on the left of the turning point. In 
this region, the general WKB solution is: 



1 



+ D2e 3 



^,»3/2 



(C.27) 



The linear solution takes the following asymptotic form as w' — )• (—1)3 x — oo: 



Cie 2" 



where 



-w 2 



2i 



3' ' 3 

Matching the linear solution to the WKB solution to the left of the turning point gives: 



-iD2 



(C.28) 

(C.29) 

(C.30) 
(C.31) 



C.3 The turning point at Z- 

We now repeat the derivation of the connection formulas at z 
region around z- where the linear approximation is valid is: 

— wM(z_,aQ,l3) 



Z-.. The general solution in the 
[CiAi{w') + C2Bi{w')] , (C.32) 
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As p is now negative, and w' ~ {—p)3w, we need the asymptotic behavior of the Airy functions for 



w 



-oo: 



1 



1 

2/^(-u;')^ ^ 



iZL 21/ /N3/2 ^ ^C_^'l3/2 



and for w' ^ oo: 



1 -2u;'3/2 

-e 3*^ 



1 2 ,3/2 

-63'" 



The WKB solution to the right of the turning point is: 

wM(z_ ,a,l3) 



1 



Die— 



2v^ ,,,3/2 



+ Doe 3 



VM,,,3/2 



where Di and D2 are related to Di and in ( P-30 ) by: 



Di = VtDi 

D2 = ^D2 
J- 1 



(C.33) 
(C.34) 

(C.35) 
(C.36) 



(C.37) 

(C.38) 
(C.39) 



Tt = e 



.^+\R'{z')\dz' 



On the left hand side of the turning point, the general WKB solution is: 



t(;M(zQ ,CK,/3) 

e 2 

1 



(C.40) 



(C.41) 



Matching the linear and WKB solutions to the left and right of the turning point as before, gives: 



Gi = -iBo ( + ^ 



Go = Bn I Yf 

z u I t 



(C.42) 
(C.43) 



Thus, using these formulas we can extend the WKB solution across the turning points, from the 
region z ^ 1 to the region z 1. 
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